Scalars convected by a 2D incompressible flow 
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1 Abstract 

We provide a test for numerical simulations, for several two dimensional 
incompressible flows, that appear to develop sharp fronts. We show that in 
order to have a front the velocity has to have uncontrolled velocity growth. 

2 Introduction 

The aim of this paper is to study the possible formation of sharp fronts in 
finite time for a scalar convected by a two dimensional divergence-free velocity 
field, with x = (x 1 ,x 2 ) G R 2 or R 2 /Z 2 } and t G [0,T) with T < oo. The 
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scalar function 8(x,t) and the velocity field u(x,t) = (ux(x,t),U2(x,t)) e R 2 
satisfy the following set of equations 

(1) 







where V^/ = (-#■ 



8x2 ' dxi 



(d t + u-V)9 

= u, 

for scalar functions f. The function ip is the stream 
function. 

There are many physical examples where the solutions satisfy the equa- 
tions above, with an extra equation or operator that relates 9 with the ve- 
locity field. Examples include; Passive scalars, Unsteady Prandtl equations, 
2D incompressible Euler equations, Boussinesq, 2D Ideal Magnetohydrody- 
namics and the Quasi-geostrophic equation. 

In the literature on numerical simulations for the 2D Ideal Magneto- 
hydrodynamics (MHD) a standard candidate for a current sheet formation 
(see Fig. 1) is when the level sets of the magnetic stream function ( repre- 
sented in (1) by 9) contain a hyperbolic saddle (an X-point configuration). 
The front is formed when the hyperbolic saddle closes, and becomes two 
Y-points configuration joined by a current sheet. (See Parker [I2|, Priest- 
Titov- Packard flT~3 1 



Friedel-Grauer-Marliani f 10 1 and Cordoba-Marliani 






t < T 



t = T 



Fig. 1. Level curves of 9. 

The same configuration was observed in numerical simulations for the 
Quasi-geostrophic equation (QG). In this case the geometry of the level sets of 
the temperature has a hyperbolic structure (See Constantin-Majda-Tabak J|], 
Okhitani-Yamada |Tl|], Cordoba || and Constantin-Nie-Schorghofer ||). The 
QG literature discusses X-points, but not Y-points. In the case of Boussi- 
nesq there is no mention, on any numerical simulation study, that a possible 
singularity is due to the closing of a hyperbolic saddle. In the work of Pumir- 
Siggia |l4j] there has been observed evidence for a formation of a front in finite 
time, across which 9 varies dramatically, on a cap of a symmetric rising bub- 
ble. E-Shu performed numerical simulations with the same initial data 
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as in which suggest that the thickness of the bubble decreases only 

exponentially. 

The equations for MHD, QG and Boussinesq are as follows 
MHD: 

(d t + u-V)6 = 
(d t + u-V)u = V ± 9-V(A9) 
u = V ± 4> 

and initial conditions 9(x, 0) = 9 and u(x, 0) = u . The V -1 ^ represents the 
magnetic field, A9 represents the current density and u = —Aip the vorticity. 
QG: 

(d t + u- V)0 = 

u = V^V where 9 = —(—A)^ip 

and initial condition 9(x, 0) = 9q. The temperature is represented by 9. 
Boussinesq: 

{d t + u-V x )9 = 
(d t + u- Va,) w = -9 X1 

u = 

Again, 9 and u are specified at time t=0. 

3 Criterion 

A singularity can be formed by collision of two particle trajectories. A trajec- 
tory X(q,t) is obtain by solving the following ordinary differential equation 

^ = u(X( g ,t),t) 
X(q,0) = q 

Therefore, 

{X(q, t) - X(p, t)) t < \X(q, t) - X(p, t) 1 1 V«| L « 
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\X(q,t) - X(p,t)\ > \X(q,0)-X(p,0)\e-SoW^ 

By this trivial argument; in order to have a collision the quantity J Q * \Vu\l°°<Is 
has to diverge. 

A classic criterion for formation of singularities in fluid flows is the the- 
orem of Beale-Kato-Majda (BKM); (see IJ), which improves the estimate 
described above, and deals with arbitrary singularities, not just collisions. 
Analogues of the BKM theorem for the above 2-dimensional equations in- 
clude the following results 

For MHD,a singularity cannot develop at a finite time T, unless we have 

f T 

/ sup x \u(x, t)\ + sup x \A x 9(x, t)\dt = oo, 
Jo 

where cu denotes the vorticity.(See Caflisch-Klapper-Steele [§].) 

For QG, a singularity cannot develop at a finite time T, unless we have 



/ sup x \ V x 9(x, t)\dt = oo, 
Jo 



(See Constantin-Majda-Tabak [|]]). 

For Boussinesq, if a singularity develops at a finite time T then 



/ sup x \u(x, t)\dt = oo and / / sup x \V x 9(x, s)\dsdt = oo. 
Jo Jo Jo 

(See E-Shu §.) 

See also Constantin-Majda-Tabak Q and Constantin-Fefferman-Majda H 
for other conditions involving direction fields, that rule out formation of sin- 
gularities in fluids. 

In the case of 2D Euler, a singularity cannot develop at a finite time. 
From the BKM viewpoint this follows from the fact that u is advected by 
the fluid, and therefore sup x \u(x,t)\ is independent of t. (See BKM [p]].) 

Instead of looking at particle trajectories we look at level curves. Because 
the scalar function 9 is convected by the flow, that implies that the level 
curves are transported by the flow. A possible singular scenario is due to 
level curves approaching each other very fast which will lead to a fast growth 
of the gradient of the scalar function. In this paper we present a variant of 
the BKM criterion for sharp front formation. We provide a test for numerical 
simulations that appear to develop sharp fronts. The BKM Theorem shows 
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that the vorticity grows large if any singularity forms; our Theorem 1 shows 
that the velocity grows large if a sharp front forms. 

The theorem we present in this paper was announced in [/J. 

4 Sharp Fronts 

The scalar function 9 is convected by the flow, therefore the level curves move 
with the flow. A sharp front forms when two of these level curves collapse on 
a single curve. We define two level curves to be two distinct time-dependent 
arcs r + (£), that move with the fluid and collapse at finite time into a 

single arc V. More precisely, suppose the arcs are given by 



r± = {fa, x 2 ) ER 2 :x 2 = f±(x, t),X\ G [a, b}} for < t < T, (2) 
with 

f ± eC\[a,b]x[0,T)) (3) 

and 

f-(x u t) < f+(xi,t) for all X!e[a,6], te[0,T). (4) 

We call the length b-a of the interval [a,b] the length of the front. The 
assumption that T±(t) move with the fluid means that 

u 2 {xi,X2,t) = -Q—(x 1 ,t)-Ui(xi,x 2 ,t) + -^-(x 1 ,t) at x 2 = f±(x u t). (5) 

This holds in particular for level curves of scalar functions g(x,t) that satisfy 
(d t + u ■ Vjb) g = 0. The collapse of T±(t) into a single curve T at time T 
means here simply that 

limt->T-(f+(xi,t) - f-(xi,t)) = for all Xie[a,b]. (6) 

and f + (xi,t) — f-(xi,t) is bounded for all X\ G [a, b], t G [0,T). 

When (g), (^), (M, (^) and (Bp hold, then we say that the fluid forms a 
sharp front at time T. 

The standard candidates for a singularity for MHD and QG are described 
by the definition given for a sharp front. We investigate the possible forma- 
tion of a sharp front. 
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The following assumption will allow us to rule out formation of sharp 
fronts. We say that the fluid has controlled velocity growth if we have 

T 

sup{\u(xi,x 2 , t)\ : xi G [a, b], f-(xi, i) < x 2 < f+(%i, t)}dt < oo. (7) 

o 

If (|7|) fails, then we say that the fluid has uncontrolled velocity growth. 

Lemma 1. Let 8 be a smooth solution of Eq.l defined for t G [0,T). Assume 
there is a sharp front at time T. Then 



b 

[f+(xi,t) - f-(x!,t)]dxi ) = if)(a,f + (a,t),t)-i/;(a,f-(a,t),t) 

+ iP(bJ-(b,t),t)-iP(bJ + (b,t),m 



Proof: Take the derivative of the stream function with respect to x\ along 
an arc T±(t) 

dj;{x 1 ,f±{x 1 ,t),t) df± 

^7 = U2{x h f±(x h t),t) - -^Mi(xi,/ ± (xi,t),t) (9) 

by combining (Q) and (|) we obtain 

ox i at 

Expression (^|) follows from integrating ( |T0D with respect to x± between a 
and b. 

Theorem 1. Let u(x,t) be a divergence-free velocity field, with controlled 
velocity growth. Then a sharp front cannot develop at time T. 

Proof: Assume there is a sharp front at time T. We define 

Mt) 

A(t)= / \J + {x x ,t) - f-{x u t)]dx x 

Ja{t) 

where 

f T 

a(t) = a + / sup{\u(xi, x 2 , s)\ : x% G [a, b], f-[xi, s) < x 2 < f+(xi, s)}ds 



and 

f-T 



b(t) = b — J sup{\u(xi, X2, s)\ : xi E [a, b], f-(xi, s) < £2 < f+{ x i> s)}ds 

There is controlled velocity growth, therefore there exists t* E [0, T) such 
that a(t) E [a, b] and b(t) E [a,b] for all t E [t*,T). 
We take the derivative of A(t) with respect to time 

dA(t) . , ... , f m d 



dt 



sup\u\ ■ S(b, t) + sup\u\ ■ S(a, t) + / — [f + (xi, t) — f-(x\, t)]dx\ 

Ja{t) Vt 



where sup\u\ = sup{\u(x\, X2, t)\ : X\ E [a, b], f-(xx, t) < x 2 < f+(xx,t)} and 
6{z,t) = f+(z,t)-f r (z,t). 

Using the definition of the stream function, the mean value theorem and 
(H), it is easy to check that > for t > t*. This contradicts (J6|) by the 
dominated convergence theorem. 

Acknowledgments 1. This work was initially supported by the American 
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